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Control of unstable waves in three-dimensional boundary layers 


Abstract 

Stability experiments are conducted in the Arizona State University Unsteady Wind 
Tunnel on a 45° swept airfoil. The pressure gradient is designed to provide purely 
crossflow-dominated transition; that is, the boundary layer is subcritical to Tollmien- 
Schlichting (T-S) disturbances. The airfoil surface is hand polished to a 0.25 fxm rms 
finish. Under these conditions, stationary crossflow disturbances grow to nonuniform 
amplitude due to submicron surface irregularities near the leading edge. Spectral de- 
compositions isolate single-mode growth rates for the fundamental and harmonic distur- 
bances. The measurements show early nonlinear growth causing amplitude saturation 
well before transition. Comparisons with nonlinear PSE calculations show excellent 
agreement in both the amplitude saturation and the disturbance mode shape. 


Background and motivation 

The present swept-wing research program experimentally investigates the fundamental 
nature of the crossflow instability which leads to transition in three-dimensional bound- 
ary layers. Gregory et al. (1955) provide the theoretical basis for the instability. It 
results in an Orr-Sommerfeld type solution that can be implemented in a variety of 
ways. See Mack (1984) for the development of the details of the instability. Reed and 
Saric (1989), Saric (1992b) and Arnal (1992) review the literature. 


Fundamentals 

In contrast to Tollmien-Schlichting (T-S) instabilities, the crossflow problem exhibits 
stationary (/ = 0) as well as traveling disturbances. The traveling waves are more 
amplified according to linear theory, however many experiments are dominated by sta- 
tionary waves. Muller and Bippes (1989), Bippes (1990, 1991), Bippes et al. (1991), 
and Deyhle and Bippes (1996) have shown that traveling waves are observed in tun- 
nels rich in unsteady freest ream disturbances, whereas stationary waves dominate in a 
low-turbulence environment. Since the flight environment is more benign than the wind 
tunnel, one expects the low-turbulence results to be more important. 
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One of the important results to come out of the DLR experiments is the set of data 
that show early saturation of the stationary disturbance amplitude and the failure of 
linear theory to predict the growth of the instability (Bippes and Nitschke-Kowsky 1990; 
Bippes et al. 1991; Deyhle et al. 1993). Dagenhart et al. (1989, 1990) and Radeztsky 
et al. (1994) observed similar behavior. Kohamaet al. (1991) showed that the stationary 
crossflow vortex controls transition by causing a high-frequency secondary instability 
resulting from the nonlinear mean-flow distortion. More recently, Deyhle and Bippes 
(1996) document the role of freestream disturbances with regard to traveling crossflow 
waves. Radeztsky et al. (1993) show that the receptivity process for stationary waves is 
strongly influenced by surface roughness, i.e., initial amplitudes. Insofar as stationary 
crossflow waves are concerned, it is clear that a successful transition prediction scheme 
must account for the initial conditions and the nonlinear growth of the disturbance. It 
is this aspect of the problem which we address in this paper. 


Objectives 

In earlier ASU experiments, Dagenhart et al. (1989, 1990) found that naturally occur- 
ring stationary crossflow waves of moderate amplitude have lower growth rates than 
predicted by linear theory. Under the same conditions, Radeztsky et al. (1993) inves- 
tigated the sensitivity to surface roughness and showed early saturation of the natural 
stationary disturbance amplitude. Later experiments by Radeztsky et al. (1994) exam- 
ined the growth of very weak crossflow waves in an attempt to close the gap between 
previous experimental results and linear theory. Even for these weak waves, linear theory 
completely failed to predict the disturbance growth. 

In all of the early experiments measuring the growth of the stationary crossflow 
wave (Dagenhart et al. 1989, 1990; Bippes et al. 1991), the initial conditions for the 
disturbance amplitude came from the unknown natural roughness of the surface. Al- 
though a dominant wavelength appears, the resulting stationary structure is nonuniform 
in span and contains many fundamental disturbance modes of unknown amplitude. Con- 
sequently, comparisons with single-mode numerical predictions are not possible. 

The same NLF(2)-0415 airfoil (Somers and Horstmann 1985) is used as in the pre- 
vious experiments. With a 45° sweep and a —4° angle of attack, the favorable pres- 
sure gradient produces considerable crossflow while suppressing T-S modes. Arrays of 
6 fim roughness elements near the leading edge produce uniform stationary disturbances 
without excessive initial amplitudes. Spectral techniques are used to identify and follow 
specific stationary modes, thus providing single- wavelength growth rates for comparison 
with theoretical calculations. 


The experiment 
Facility and model 

The Unsteady Wind Tunnel at Arizona State University is a low-speed, low-turbulence, 
closed-circuit facility used to study the stability and transition of laminar boundary 



layers (Saric 1992a). The airfoil model is mounted vertically in the 1.4 m x 1.4 m x 5 m 
test section. The aluminum surface is hand polished to a 0.25 fm x rms finish so that 
even micron-sized roughness elements are well above the background roughness level. 


Measurement techniques 

Two standard hot-wire scanning methods are used to investigate the stationary crossflow 
waves. These are described below. 


Wall-normal scans 

Wall-normal boundary-layer scans provide a detailed map of the stationary structure. 
These maps are constructed by taking a spanwise series of mean-flow boundary-layer 
profiles at constant x/c. Once the scans are aligned with the airfoil surface, disturbance 
profiles are generated from which a stationary crossflow mode shape is computed. The 
disturbance growth is calculated by tracking the size of the mode shape at various chord 
positions. Three different measures of the mode shape size are used to characterize the 
disturbance amplitude: the maximum of the mode, the integral of |it ; | with respect to 
y, and the integral of \u'\ 2 with respect to y. When the disturbance growth is cast in 
terms of the amplification factor N, all three measures collapse onto a single iV-factor 
curve. These results, however, cannot (in general) be quantitatively compared with 
single-wavelength linear predictions since all amplified stationary modes are lumped 
into a single mode shape. 

Single-wavelength information can be extracted from a boundary-layer map by taking 
a spanwise slice across the profiles at a constant height above the airfoil surface. The 
resulting velocity vs. span trace can be decomposed using spectral techniques to resolve 
the wavenumber content. Reasonable resolution in the wavenumber domain, however, 
requires a large spanwise extent of the measurement region. Consequently, this technique 
can quickly become prohibitively time consuming. 


Spanwise scans 

Individual- wavelength growth rates are obtained by restricting hot-wire measurements to 
a single spanwise scan at a constant height above the airfoil surface. With this technique, 
data are acquired along the entire span of the measurement region (240 mm) at much 
higher (spanwise) resolution than the wall-normal scans. Since data are collected at only 
one height in the boundary layer, these scans progress very quickly (typically 75 minutes 
per scan compared to nearly 45 hours for a full set of wall-normal scans with the same 
spanwise extent and resolution). The disadvantage of this technique is that the details 
of the stationary structure are not captured. For this reason, a disturbance mode shape 
(produced with a small set of wall-normal scans) is used to guide the spanwise scans. 



Results 


Basic state 

The basic state is documented with pressure measurements and mean-flow boundary- 
layer profiles. Figure 1 shows the measured and theoretical C p distribution over the 
airfoil upper surface. The experimental data are the average readings from two sets of 
pressure ports in the airfoil surface. The theoretical curve is the inviscid solution from 
the NASA Langley code MCARF. The agreement, especially in terms of the pressure 
gradient, is good over the entire measurement region. 


Boundary-layer maps 

Figure 2 is a contour plot of the streamwise velocity u/U e in the (y, z) plane. The flow 
is toward the reader, and the stationary vortices are turning in the right-handed sense. 
These data are acquired at x/c = 0.45, Re c = 3.0 x 10 6 , with no artificial roughness on 
the airfoil. The nonuniformity of the naturally occurring stationary waves is caused by 
submicron surface irregularities near the leading edge. Figure 2 displays a strong feature 
at a 12 mm spacing, which is approximately the most amplified stationary wavelength. 
At the same time, the richness in the spectral content is evident. This is typical of all 
of the earlier data (Bippes and Nitschke-Kowsky 1990; Bippes et al. 1991; Deyhle et al. 
1993; Dagenhart et al. 1989, 1990), and indicates both nonlinear behavior and multiple 
modes. Thus, even a nonlinear calculation that included only a single spanwise mode 
would be inappropriate to characterize the disturbance motion. 

In order to generate spanwise-uniform stationary crossflow waves, initial conditions 
are controlled by applying a full-span array of roughness elements at x/c = 0.023 (near 
the neutral point of the instability) following Radeztsky et al. (1993). The roughness 
height is k — 6 /im, and the 3.7 mm diameter elements are spaced 12 mm apart in span, 
corresponding to the most amplified wavelength according to linear theory. Figure 3 
shows the streamwise velocity contours at x/c = 0.45 and Re c = 2.4 x 10 6 with this 
roughness distribution. Under these conditions, the roughness Reynolds number, Rek , 
is 0.1. The uniformity of the fundamental 12 mm mode is striking, and allows for 
meaningful comparisons with single-mode theoretical predictions. A single vortex is 
isolated in Figure 4 and plotted on a 1:1 scale. The crossflow vortex produces regions of 
upwelling and downwelling which transport low- and high-momentum fluid, respectively. 
The symmetry of the co-rotating vortex distorts this momentum transfer giving an 
apparent rollover of low-momentum fluid that appears above high-momentum fluid. 

Figure 5 shows the 100 boundary-layer profiles from which Figure 3 is generated. The 
profiles are obtained at 1 mm intervals in the swept span direction. The dots indicate 
the spanwise average of the profiles, which accounts for basic state plus the mean-flow 
distortion [(0,0) mode]. It should be emphasized that these are mean profiles and not 
an unsteady oscillation in the boundary layer. One can clearly see how the stationary 
vortex structure has distorted the mean flow, resulting in accelerated, decelerated, and 
doubly inflected profiles existing millimeters apart. The nonlinearities are indicated by 



the distortion of the averaged profile in the vicinity of y « 2.8 mm. This distortion of the 
basic state leads to the secondary instability which controls the transition to turbulence 
(Kohama et al. 1991). 

Disturbance profiles are generated by subtracting the basic state plus mean-flow 
distortion (i.e., the span wise average profile) from the individual boundary-layer profiles 
(Figure 6). From these data, the crossflow mode shape is generated by computing the 
spanwise rms of the disturbance profiles (Figure 7). This mode shape contains the 
fundamental disturbance and all amplified harmonics. 

The total disturbance amplitude is computed using three measures of the mode shape 
size as outlined above. Growth rates are obtained by repeating this procedure at several 
chord positions. Figure 8 shows the disturbance amplitude distribution for 6 fim rough- 
ness with 12 mm spacing at Re c = 2.4 x 10 6 . The dashed lines represent the absolute 
size of the mode shape as computed by each of the three measures. The solid lines 
show the corresponding amplification factor N . The ability of the iV-factor to collapse 
the data onto a single curve is typical. The nonlinear saturation is clearly evident, and 
occurs well before the transition location of (x/c ) tr = 0.52. It is worth emphasizing that 
these measures record the total disturbance amplitude since all amplified wavelengths 
are present in the mode shape. 

Figure 9 compares the experimental 7V-factor (as computed from the maximum of 
the mode shape profiles) with various theoretical predictions of Haynes and Reed (1996) . 
The nonlinear parabolized stability equations (NPSE) results are computed using ini- 
tial amplitudes provided by the experiment. The agreement is excellent, especially in 
predicting the amplitude saturation. (At this time, the Haynes and Reed formulation 
does not contain curvature, which is known to be stabilizing and may account for the 
small differences in the disturbance growth.) In contrast, the Orr-Sommerfeld and linear 
PSE results fail to predict the details of the disturbance growth. The early qualitative 
agreement with the linear PSE results indicates that the nonlinear effects are initially 
weak up to x/c = 0.20, at which time the growth rates depart from linear behavior. It 
is at this location that the spanwise average of the boundary-layer profiles first begins 
to exhibit the distortion shown in Figure 5. 


Wavelength separation 

As pointed out previously, crossflow amplitudes computed from mode-shape profiles 
contain all amplified stationary modes. To decompose the wavenumber content of the 
disturbance, the spanwise scan technique discussed above is used. Figure 10 shows the 
spanwise velocity profile at x/c = 0.45 for Re c = 2.4 x IQ 6 . The roughness configuration 
is k = 6 pm with 12 mm spanwise spacing. The data are acquired at y = 0.9 mm 
above the airfoil surface, corresponding to the maximum of the mode-shape profile. The 
strong distortion of the boundary-layer is evident, as is the spanwise uniformity of the 
stationary structure. 

The FFT-based power spectral density (PSD) is shown in Figure 11. The dominance 
of the fundamental 12 mm mode is clear. However, the super harmonics at A s = 6 mm 
and 4 mm are also amplified. (The amplitude of the 4 mm mode is too small to appear 



on the scale of Figure 11.) This wave doubling was observed by Saric and Yeates (1985) 
and predicted by Reed (1988). Disturbance amplitudes are computed by integrating 
the peaks of the PSD. When repeated at several chord positions, the growth rates of 
Figure 12 are obtained. The dashed lines represent the rms amplitude of the individual 
modes, while the amplification factor N is plotted with solid lines. At each chord 
location, the spanwise scan is taken at the boundary-layer height corresponding to the 
maximum of the mode-shape profile. This permits us to compare the amplitude of the 
individual modes with that of the total disturbance as computed from the maximum of 
the mode shape (indicated on Figure 12 by the lines marked “Total”). It is a remarkable 
verification of the two different measurement techniques that the amplitudes of the total 
and fundamental (X s = 12 mm) disturbances agree for x/c < 0.30. For x/c > 0.30, the 
amplitude of the fundamental diverges from the total disturbance, and the harmonic 
at X s = 6 mm becomes measurable. In the region of strong nonlinear interaction and 
amplitude saturation (x/c > 0.35), the 6 mm component shows significant energy and 
the 4 mm mode becomes unstable. 


Conclusions 

Stationary crossflow waves are investigated on a swept airfoil within a low-disturbance 
environment. It is shown that even though the surface is polished to 0.25 fim rms 
roughness, stationary crossflow waves still dominate the transition process. Because the 
initial roughness is nonuniform, the resulting disturbance motion is nonuniform and is 
complicated by the presence of many different modes. Because of this, comparisons with 
different theoretical and computational schemes are unnecessarily rendered much more 
difficult. 

Systematic introduction of equally spaced 6 /im roughness elements (Rek = 0.1) is 
shown to produce an ideal fundamental mode at that spacing. When the roughness 
elements are placed at the most amplified linear-mode wavelength of X s = 12 mm, 
disturbance growth, departure from linearity, and saturation amplitude axe documented. 
Evidence of growth at the 6 mm and 4 mm harmonics is shown. There is no evidence 
of wavelength doubling (i.e., components at 24 mm or 36 mm). Comparisons with 
computations using nonlinear parabolized stability equations are very good. It is possible 
to conclude that the introduction of systematic weak roughness provides the necessary 
data base for comparisons with theory and computations. 

These data continue to illustrate the extreme sensitivity to leading-edge roughness. 
Strong nonlinear distortion of the mean flow is observed, as is the nonlinear saturation 
of the disturbance amplitude. This saturation occurs well before the transition to tur- 
bulence. Linear theory fails to capture these details of the disturbance growth. On the 
other hand, NPSE calculations agree remarkably well with the experimental data. 
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Figure 1: Experimental and theoretical swept 
C p distribution at Re c = 2.4 x 10 6 . The the- 
ory is computed with the NASA Langley code 
MCARF. 
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Figure 3: Streamwise velocity contours at 

x/c = 0.45, Re c — 2.4 x 10 6 . A full-span array 
of 6 (im roughness with 12 mm spacing is at 
x/c = 0.023. 
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Figure 2: Streamwise velocity contours at 

x/c = 0.45, Re c == 3.0 x 10 6 . No artificial 
roughness. 


Figure 4: Single stationary crossflow vortex 
isolated from Figure 3 and plotted on a 1:1 
scale. 
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Figure 5: Spanwise array of boundary-layer 
profiles at x/c = 0.45, Re c = 2.4 x 10 6 . A 
full-span array of 6 fim roughness with 12 mm 
spacing is at x/c = 0.023. The dots indicate 
the mean of the profiles. 
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Figure 6: Spanwise array of disturbance pro- 
files at x/c = 0.45, Re c = 2.4 x 10 6 . A full- 
span array of 6 /m l roughness with 12 mm 
spacing is at x/c = 0.023. 


Figure 7: Stationary crossflow mode shape at 
x/c = 0.45, Re c — 2.4 x 10 6 . A full-span array 
of 6 fJLm roughness with 12 mm spacing is at 
x/c = 0.023. 



Figure 8: Disturbance amplitude distribution 
at Re c = 2.4 x 10 6 . A full-span array of 6 /xm 
roughness with 12 mm spacing is at x/c = 
0.023. The reference point for the N-factor 
calculations is x/c = 0.10. 










Figure 9: Measured and theoretical N-i actor 
for Re c == 2.4 x 10 6 . A full-span array of 6 fim 
roughness with 12 mm spacing is at x/c = 
0.023. The reference point for the iV-factor 
calculations is x/c ■■== 0.10. 



Figure 11: FFT power spectrum of span wise 
hot-wire scan at x/c = 0.45, #e c == 2.4 x 10 6 , 
y = 0.9 nun. A full-span array of 6 rough- 
ness with 12 mm spacing is at x/c = 0.023. 








